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ABSTRACT: The most general description of the classical world is in terms of local 
densities (such as number, momentum, energy), and these typically evolve according to 
evolution equations of hydrodynamic form. To explain the emergent classicality of these 
variables from an underlying quantum theory, it is therefore necessary to show, firstly, 
that these variables exhibit negligible interference, and secondly, that the probabilities for 
histories of them are peaked around hydrodynamic evolution. The implementation of this 
programme in the context of the decoherent histories approach to quantum theory is de- 
scribed. It is argued that, for a system of weakly interacting particles, the eigenstates of 
local densities (averaged over a sufficiently large volume) remain approximate eigenstates 
under time evolution. This is a consequence of their close connection with the correspond- 
ing exactly conserved (and so exactly decoherent) quantities. The subsequent derivation 
of hydrodynamic equations from decoherent histories is discussed. 



If the universe is described at the most fundamental level by quantum theory, why 
is it so very nearly classical? There are very many aspects to the issue of emergent 
classicality (see, for example, Ref.[l] for an overview), but crucial to most of them is 
the demonstration that certain types of quantum states of the system in question exhibit 
negligible interference. Initial superpositions of such states may therefore be effectively 
replaced by statistical mixtures. This, loosely speaking, is decoherence, and has principally 
been demonstrated for the situation in which there is a distinguished system, such as a 
particle, coupled to its surrounding environment [2,3]. 

Most generally, decoherence typically comes about when the variables describing the 
entire system of interest naturally separate into "slow" and "fast", whether or not this 
separation corresponds to, respectively, system and environment t . If the system consists 
of a large collection of interacting identical particles, as in a fluid for example, the natural 
set of slow variables are the local densities: energy, momentum, number, charge etc.. These 
variables, in fact, are also the variables which provide the most complete description of 
the classical state of a fluid at a macroscopic level. 

The most general demonstration of emergent classicality therefore consists of showing 
that, for a large collection of interacting particles described microscopially by quantum 
theory, the local densities become effectively classical. Although one might argue that 
the system-environment mechanism might play a role, since the collection of particles are 
coupled to each other, decoherence comes about in these situations for a different reason: 
it is because the local densities are almost conserved if averaged over a sufficiently large 
volume [5]. Hence, the approximate non-interference of local densities is due to the fact 
that they are close to a set of exactly conserved quantities, and exactly conserved quantities 

t See Ref.[4] for a discussion of the conditions under which the total Hilbert space may be 
written as a tensor product of system and environment Hilbert spaces 
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obey superselection rules. 

Intuitively appealing though this argument is, it is clearly a quantitative issue. The 
object of this letter is to show that, under certain reasonable conditions, local densities 
averaged over a sufficiently large volume are indeed approximately decoherent as a result 
of their close connection to exact conservation. 

We will approach the question using the decoherent histories approach to quantum 
theory [5,6,7,8], which has proved particularly useful for discussing emergent classicality 
in a variety of contexts. The central object of interest is the decoherence functional, 



The histories are characterized by the initial state and by the strings of projection op- 
erators P a at times t\ to t n (and a denotes the string of alternatives a\ • • • a n ). Intuitively, 
the decoherence functional is a measure of the interference between pairs of histories a, 
a'. When it is zero for we say that the histories are decoherent and probabilities 

p(a) = D(a,a) obeying the usual probability sum rules may be assigned to them. One 
can then ask whether these probabilities are strongly peaked about trajectories obeying 
classical equations of motion. For the local densities, these equations will be hydrodynamic 
equations, and these and closely related aspects of emergent classicality have been pursued 
at greater length elsewhere Refs. [4,9,10]. 

We consider the class of systems which are described at the microscopic level by a 
Hamiltonian of the form 





(1) 
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For definiteness, we will concentrate on the case of a weakly interacting dilute gas, making 
brief reference to a one-dimensional chain of oscillators, but it will be clear that the physical 
ideas are reasonably general. The local densities of interest are the number density n(x), 
the momentum density g(x) and the energy density ft(x), defined by, 

n(x) = £ *(x- qj ) (3) 
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gW = Pi ^-qj) ( 4 ) 

j \ e>i J 

(suitably ordered, in the quantum case). We are interested in local densities smeared over 
a volume V. The effect of this is to replace the delta functions with a window function, 
denoted 5y, which is zero outside V and 1 inside. It is also useful to work with the Fourier 
transforms of the local densities, denoted n(k), g(k), h(k). So, for example, the local 
number density at wavelength k is 

n(k) = e' kqb (6) 
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Exact conservation is obtained in the limit k = |k| — > 0, or V — > oo in (3)-(5). 

We would like to compute the decoherence functional for histories consisting of projec- 
tions onto the operators (3)-(5). (The construction of the projectors is described in more 
detail in Ref.[10]). In the case of exact conservation, k = 0, we have exact decoherence 
simply because the projectors in Eq.(l) all commute with H and with each other [11]. 
Our main task is therefore to show that as k increases from zero there is still a non-trivial 
regime in which decoherence is approximately maintained. 

A significant result in this direction has been established already by Calzetta and Hu 
for the case of local temperature T(x) obeying the diffusion equation [12]. They took their 
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initial state to be close to the equilibrium state, and worked backwards from the diffusion 
equation plus fluctuations to deduce the influence functional it must have arisen from, from 
which the degree of decoherence could be deduced. Here, by contrast, initial macroscopic 
superposition states are considered. A more detailed comparison of these two approaches 
is certainly of interest. 

We begin by rewriting the exact conservation case in a simple way that makes its 
generalization to locally conserved quantities more apparent. Suppose the histories are 
projections onto some conserved quantity, Q. Let the initial state be a superposition of 
eigenstates of Q, 

|*> = ^ (|a) + |6» (7) 

where (a\b) = and 

Q\a) = a\a), Q\b) = b\b) (8) 

Since the P a 's are projections onto Q, P a either annihilates or preserves \a) and \b). Take 
the case of a history with just two moments of time (the generalization to more times is 
trivial). The only non-zero off-diagonal terms of the decoherence functional are of the form 

D{a,g[) = ^Tr (p a2 e-^ m \a){b\e^ Hi ) 

= l -Tt{P a2 \a t ){b t \) (9) 

But Q is conserved, hence [P a2 , H] = and 

Pa 2 \a t ) = Pa 2 e-i Ht \a) 

= e -i Ht P a2 \a) = \a t ) (10) 

(or equals zero if does not correspond to a). It follows that 

D(a, a') = ^Tr (P a2 \a t ) (b t \) 

= (b t \a t ) = (b\a) = (11) 
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and therefore we have decoherence. 

Now suppose that the operator Q is one of the local densities (3)-(5), so is no longer 
exactly conserved. The steps up to Eq.(9) still hold. But to go further, we need to 
know how the eigenstates of the local densities behave under time evolution. A reasonable 
supposition, which will be justified, is the following. Let us suppose that under time 
evolution, the eigenstates of Q remain approximate eigenstates. That is, we initially have 
(7), but under evolution to time t, 

QWt) « (Q)\at) (12) 

or, more precisely, 

^ « 1 (13) 

i.e., the state remains strongly peaked in the variable Q under time evolution. The states 
are then approximate eigenstates of the projectors, so that in place of Eq.(lO), we have the 
approximate result, P a2 \at) ~ \o>t) ( or equals zero) as long as the width of the projection 
is much greater than the uncertainty (AQ) 2 . Hence Eq.(ll) follows approximately, and 
we get approximate decoherence to the extent that the approximation (13) holds. 

The key point is therefore the following: approximate decoherence is assured for histo- 
ries of operators Q whose eigenstates have the property that they remain strongly peaked 
in Q under time evolution, as characterized by (13). To demonstrate decoherence of the 
local densities, therefore, we need only find their eigenstates, and show that they satisfy 
the localization property (13) under time evolution. (Note, incidently, that the above argu- 
ment actually assures decoherence of any variables Q satisfying the localization property. 
The particular significance of the local densities is that they are continuous functions of the 
coarse graining scale k, so are guaranteed to satisfy the requisite property if k is sufficiently 
close to zero.) 
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Since the three operators (3)-(5) do not commute, exact simultaneous eigenstates 
cannot be found. However, there are approximate simultanous eigenstates. For weak 
interactions, they are products of N identical terms, 

I*) = W) ® |V>) ® • • • ® W) (14) 

and are approximate eigenstates of all three operators for large N. The proof of this 
statement involves considering, for the local number density for example, the object 
(An(x)) 2 /(n(x)) 2 , and showing that it goes like 1/N for large N (see Ref.[10], for ex- 
ample). It is essentially the central limit theorem (see also Ref.[13]). For the number 
and momentum density it relies on the fact that they are sums of identical one-particle 
operators. For the local energy density, it additionally requires the smearing volume to 
be sufficiently large, compared to some lengthscale indicated by the interactions. Some 
tuning of the state can be carried out to ensure that (14) is an optimal approximate 
eigenstate of all the local densities but this will not be done here. (Also, the passage to 
exact eigenstates of n(k), <7(k), h(k) as k — > can be seen explicitly if the one-particle 
states IV') are taken to be one-particle momentum eigenstates). 

The question is now what happens to the eigenstates (14) of the local densities under 
time evolution by the Hamiltonian (2). Consider first the trivial but enlightening case 
in which there no interactions. In this case, the time evolved eigenstates \at) remain of 
the product form (14), so they are still approximate eigenstates of the local densities (but 
with a time-evolved eigenvalue) for the same reasons as above. Hence there is approximate 
decoherence. 

Decoherence in the non-interacting case comes about for two reasons. First, it is due 
to the fact that a state of the form (14) will remain strongly peaked about the average 
values of the local densities, n(x), g(~x), h(x.) under time evolution, and thus the state is 
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essentially undisturbed by the projectors (as long as their widths are sufficiently large). 
The strong peaking follows from the assumption of large N and from the fact that the 
local density operators are sums of identical one particle operators. Secondly, it is due 
to the almost trivial fact that the orthogonality of the two elements of the initial state is 
preserved by unitary evolution. 

This second fact is important because the first one is not always sufficient to guarantee 
decoherence. Although the state remains strongly peaked about the average values of the 
local densities, these average values do not necessarily obey deterministic equations. In 
the case of histories characterized by number density only, for example, (n(x)) at time t 
is not uniquely determined by (n(x)) at the initial time (in the state (14)). That is, in 
Eq.(9), | at) and \bt) may in fact be peaked about the same value of number density, even 
though the initial values are different. The decoherence is therefore not in fact due to an 
approximate determinism (such as that used in the phase space histories of Omnes [7]). It 
is necessary only that the evolved states are essentially undisturbed by the projectors, and 
therefore that the two orthogonal components of the initial state are eventually overlapped 
at the final time, as in Eq.(ll), to give zero. 

The next and most important task is to show that the above story is in fact still true, 
with qualifications, in the presence of interactions. The complete description of N interact- 
ing particles is generally extremely involved, but we can make some progress by restricting 
attention to a sufficiently dilute gas of weakly interacting components, and then making 
two assumptions which are standard in kinetic theory and non-equilibrium statistical me- 
chanics [14] . It is notationally convenient in what follows to work with a Wigner function, 
rather than quantum state. Hence associated with the full iV-particle wave function is 
an iV-particle Wigner function Wjv(Pl) qi, • • 'Piv, qjv)- For a dilute, weakly interacting 
gas, it is reasonable to assume that three-particle correlations are negligible. This is our 
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first assumption. It means that all the physics is contained in the one and two-particle 
reduced Wigner functions, Wi(pi,qi) and W^(pi, qi, P2, q2)- All higher order reduced 
Wigner functions will reduce to products of these. 

We again take as our initial state the product state (14) (which is still an approximate 
eigenstate in the interacting case), and let it evolve, so correlations will develop. The degree 
to which the particles become correlated is contained in the two-particle distribution W2 of 
the evolved eigenstate. On general grounds, we expect that the inter-particle correlations 
will only be important on some length scale L, and beyond that length scale, they will be 
uncorrelated. That is, we will assume that 

W2(Pl,qi,P2,q2) ~ Wi(pi,qi)Wi(p2,q2) (15) 

for |q2 — qi| > L, and otherwise W2 will have a form indicating non-trivial correlations. 
This is our second assumption. It is physically reasonable, and it is in fact a key assumption 
in the derivation of the Boltzmann equation [14] . 

Note that the assumption (15) would not necessarily be appropriate for all possible 
initial quantum states. One could construct initial quantum states which would possess 
or develop non-trivial long-range correlations, for which this assumption may never hold. 
However, as seen in the non-interacting case, to demonstrate decoherence we only need to 
consider the time evolution of the special class of initial states which are eigenstates of the 
local densities. These particular initial states, which are of the form (14), approximately, do 
not have long-range correlations. It is therefore very plausible, at least for a dilute, weakly 
interacting gas, that they will develop only limited correlations under time evolution and 
the assumption (15) will hold. 

Given the above assumptions, it is now reasonably straighforward to argue that the 
state is still strongly peaked about the average values of the local densities, as long as 
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V » L 3 . For example, for the number density, we have 

(n(x)> = ^(Sviqj - x)> = N f d\ p(q) (16) 

where p(q) is the one-particle probability distribution of q (obtained by integrating the 
one-particle Wigner function over p). Similarly, 

(n 2 (x)) = J2 (MQj - *)Mq« - x)) 

= N(S V ) + (N 2 - N)(S v ( qi - x)dy(q 2 ~ x)) (17) 

where we have used 8y = 8y, and also an assumption of identical particles to reduce the 
sum over j, £ to particles labeled 1 and 2. We now have 

(An(x)) 2 =<n 2 (x)> - (n(x)) 2 

=N 2 ((^(qi - x)Mq 2 " x)> " (V) 2 ) 

+ N ((Sy) - (5y( qi - x)Mq 2 " x))) (18) 

If there is no correlation at all between the particles, the coefficient of iV 2 would vanish, 
so (An(x)) 2 /(n(x)) 2 would go like 1/N, which goes to zero as iV — > oo. This is the 
standard central limit theorem result indicated earlier for the non-interacting case. With 
interactions, the coefficient of iV 2 is no longer zero. We now need to show, therefore, 
that this term is still sufficiently small for (An(x)) 2 /(n(x)) 2 to remain small as iV — > oo. 
Introducing the two-particle distribution p(qi,q2) (obtained by integrating pi,p 2 out of 
W2), it is readily shown that the leftover terms as — > 00 are 

(An(x)) 2 = f v rf 3 qi J v d\ 2 (p(qi, qg) ~ P(qi)p(q2)) (1Q) 

This is clearly zero if there are no correlations. In the interacting case we use the assump- 
tion (15), which implies that 

p(qi»q2) ~p(qi)p(q2) (20) 
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for |qi — q2 1 > L, and otherwise non-trivial correlations exist. Hence the integral in the 
numerator takes contributions only from the region — c&\ < L. 

To see that (19) is small, note that in the numerator, the integral is over a volume V 2 
in the six-dimensional two particle configuration space. If V « L?, the factorization of 
p(qi) O2) f° r I Ol ~~ 02 1 > L makes no difference, since qi and c\2 can never be far enough 
apart in the integrand (assuming V is regular in shape). However, if V » L 3 , the V 2 - 
sized integration region is substantially reduced in size to Fxil On dimensional grounds 
the numerator is thereofre proportional to a number of order FL 3 , and the denominator 
to V 2 (perhaps with other factors common to both). This means that 

(An(x)) 2 Z/> 

^F~v (21) 

This order of magnitude estimate becomes exact if we assume that the probabilities are 
constant in the region of non-trivial correlation (another common assumption of kinetic 
theory [14]). Hence the state will be strongly peaked about the average of n(x) if V » L 3 . 

In the one-dimensional oscillator chain model considered in Ref.[15], the uncertainty 
in n(k) (the one-dimensional version of Eq.(6)), can be computed explicitly in the special 
case of a Gaussian state. It is, 

N N 

(An(k)) 2 = J2J2^ kqj ^ e ~ ikq ^ (e PfTfe '' %) - l) (22) 
j=ll=l 

where cr(qj,q£) = (qjqg) — (qj)(qg) measures the degree of correlation between different 
particles in the chain. As k increases from zero, the leading order terms in (22) are of the 
form, (An(k)) 2 = k 2 (AX) 2 , where X = (the centre of mass coordinate), and since 

(n(k)) ~ N, we have 

(An(k)) 2 k 2 (AX) 2 



\(n(k)}\ 2 N 2 
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(23) 



This will be very small as long as k^ 1 is much larger than the lengthscale of a single 
particle. (An(k)) starts to grow very rapidly with k, and (23) is no longer valid, when 
k~ l becomes less than the correlation length indicated by cr(qj,qi)- Hence the state is 
strongly peaked about the mean as long as the coarse graining lengthscale remains 
much greater than the correlation length of the time-evolved local density eigenstates. This 
correlation length is considered in Ref.[15] and found to be generally very small compared 
to the system size. A simple field theory model is also considered in Ref.[15], confirming 
many of the expected features outlined here. 

It is possible to see on physical grounds why one expects a result of the form (21) to hold 
quite generally. In the non-interacting case we used the central limit theorem result that 
(An) 2 / (n) 2 goes like 1/N. In the interacting case, the state is no longer of the product form 
(14), but an analagous result still holds. The point is that the correlations that develop 
extend only over a (typically small) volume of size L 3 , so the system breaks up into a large 
number of essentially identical uncorrelated regions of this size. Therefore each smearing 
volume V, if much greater than L 3 , contains of order V/L^ identical uncorrelated regions 
each of which contribute equally to the local density averaged over V. Loosely speaking, 
a central limit theorem-type result again applies, not to the N uncorrelated particles in 
the same state, but to the V/L^ uncorrelated regions. So 1/N is replaced by L^/V in the 
central limit theorem, and hence the above result. 

Similar results hold for the local momentum and energy density. We have therefore 
demonstrated the desired result: eigenstates of the coarse-grained local densities remain 
approximate eigenstates under time evolution as long as the smearing volume is much 
greater than the correlation volume of these states. Decoherence of these variables then 
follows. 

We now briefly consider the probabilities for histories. They are strongly peaked at 
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each moment of time about the average values, (n(x, £)), (<7(x, £)), (h(x,t)), averaged in a 
local density eigenstate. The hydrodynamic equations (or even a closed set of equations) do 
not necessarily follow, however, since these require a local equilibrium initial state [16,14], 
but we outline how this might come about. 

The averages of the local densities depend only on the one-particle Wigner function 
(with a small correction depending on the two particle function in the case of local energy 
density), hence the evolution equation of the local densities can be determined by obtaining 
an evolution equation for the one-particle Wigner function W\. With a local density 
eigenstate as initial state, and with the two assumptions utilized above, we expect this 
evolution equation is the Boltzmann equation. We therefore evolve W\, subject to the 
initial condition that it be equal to the one-particle Wigner function of a local density 
eigenstate. This initial Wigner function is not of local equilibrium form, but it is reasonable 
to expect that it will rapidly approach local equilibrium form under evolution according 
to the Boltzmann equation, and thereafter retain that form. (The temperature, chemical 
potential etc., of the local equilibrium state will be determined by the average values 
(/i(x, t)) etc.). Hence, except for a short initial period during which the initial state 
settles down to local equilibrium form, the probabilities for histories will be peaked about 
hydrodynamic equations. 

Some of these features can be seen in some detail in Ref.[10], where the emergence of 
the diffusion equation was considered. The system studied was a collection of N foreign 
non-interacting particles in a background fluid. Decoherence was therefore provided largely 
by the fluid in this case, rather than by conservation, but the interest of the model is that 
it gives an explicit picture of the emergence of a hydrodynamic equation. Each foreign 
particle behaves like a quantum Brownian particle, whose evolution equation is well-known. 
An initial state of the form (14) for the N foreign particles evolves into a mixed state of 
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the form of an AT-fold product, 



P = PI® Pi---® Pi (24) 

where each 1-particle density operator p\ describes quantum Brownian motion. From the 
Wigner function of p\ it is readily shown that the 1-particle position distribution obeys 
the diffusion equation at long times, from which it readily follows that the iV-particle 
number density n(x) also obeys the diffusion equation. This model is therefore a kind of 
"half-way house" between the decoherence-through-environment and decoherence-through- 
conservation models, but it helps to complete the general picture. 

Summarizing, the final picture we have is therefore as follows. An initial state consisting 
of a superposition of local density eigenstates may be treated as a mixture of the same 
states, since they are decoherent. Each state separately will give probabilities peaked 
about hydrodynamic equations, with particular values of initial values of phenomenological 
parameters such as temperature etc., and these will be different for each element of the 
mixture. We therefore have a statistical mixture of trajectories, each evolving according 
to hydrodynamic equations but with different phenomenological parameters, i.e., to very 
different macroscopic states. More details of this work may be found elsewhere [15,17]. 
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